Genetic investigations on backfat thickness and body condition score in German Holstein cattle.

Up to now, little has been known about backfat thickness (BFT) in dairy cattle. The objective of this study was to investigate the lactation curve and genetic parameters for BFT as well as its relationship with body condition score (BCS) and milk yield (MKG). For this purpose, a dataset was analysed including phenotypic observations of 1929 German Holstein cows for BFT, BCS and MKG recorded on a single research dairy farm between September 2005 and December 2022. Additionally, pedigree and genomic information was available. Lactation curves were predicted and genetic parameters were estimated for all traits in first to third lactation using univariate random regression models. For BCS, lactation curves had nadirs at 94 DIM, 101 DIM and 107 DIM in first, second and third lactation. By contrast, trajectories of BFT showed lowest values later in lactation at 129 DIM, 117 DIM and 120 DIM in lactation numbers 1 to 3, respectively. Although lactation curves of BCS and BFT had similar shapes, the traits showed distinct sequence of curves for lactation number 2 and 3. Cows in third lactation had highest BCS, whereas highest BFT values were found for second parity animals. Average heritabilities were 0.315 ± 0.052, 0.297 ± 0.048 and 0.332 ± 0.061 for BCS in lactation number 1 to 3, respectively. Compared to that, BFT had considerably higher heritability in all lactation numbers with estimates ranging between 0.357 ± 0.028 and 0.424 ± 0.034. Pearson correlation coefficients between estimated breeding values for the 3 traits were negative between MKG with both BCS (r = -0.245 to -0.322) and BFT (r = -0.163 to -0.301). Correlation between traits BCS and BFT was positive and consistently high (r = 0.719 to 0.738). Overall, the results of this study suggest that BFT and BCS show genetic differences in dairy cattle, which might be due to differences in depletion and accumulation of body reserves measured by BFT and BCS. Therefore, routine recording of BFT on practical dairy farms could provide valuable information beyond BCS measurements and might be useful, for example, to better assess the nutritional status of cows.

Partitioning variance in a signaling trade-off under sexual selection reveals among-individual covariance in trait allocation.

Understanding the evolution of traits subject to trade-offs is challenging because phenotypes can (co)vary at both the among- and within-individual levels. Among-individual covariation indicates consistent, possibly genetic, differences in how individuals resolve the trade-off, while within-individual covariation indicates trait plasticity. There is also the potential for consistent among-individual differences in behavioral plasticity, although this has rarely been investigated. We studied the sources of (co)variance in two characteristics of an acoustic advertisement signal that trade off with one another and are under sexual selection in the gray treefrog, Hyla chrysoscelis: call duration and call rate. We recorded males on multiple nights calling spontaneously and in response to playbacks simulating different competition levels. Call duration, call rate, and their product, call effort, were all repeatable both within and across social contexts. Call duration and call rate covaried negatively, and the largest covariance was at the among-individual level. There was extensive plasticity in calling with changes in social competition, and we found some evidence for among-individual variance in call rate plasticity. The significant negative among-individual covariance in trait values is perpendicular to the primary direction of sexual selection in this species, indicating potential limits on the response to selection.

Short communication: Correlation of methane production, intensity, and yield with residual feed intake throughout lactation in Holstein cows.

The environmental impact of dairy production can be reduced in several ways, including increasing feed efficiency and reducing methane (CH4) emissions. There is no consensus on their relationship. This study aimed at estimating the correlations between residual feed intake (RFI) and CH4 emissions expressed in g/d methane production (MeP), g/kg of fat- and protein-corrected milk methane intensity (MeI), or g/kg of DM intake methane yield (MeY) throughout lactation. We collected CH4 data using GreenFeed devices from 107 Holstein cows, as well as production and intake phenotypes. RFI was predicted from DM intake, fat- and protein-corrected milk, BW, and body condition score. Five-trait random regression models were used to estimate the individual variance components of the CH4 and production traits, which were used to calculate the correlations between RFI and CH4 traits throughout lactation. We found positive correlations of RFI with MeP and MeI ranging from 0.05 to 0.47 throughout the lactation. Correlations between RFI and MeY are low and vary from positive to negative, ranging from -0.18 to 0.17. Both MeP and MeI are favorably correlated with RFI, as is MeY during the first half of lactation. These correlations are mostly favorable for genetic selection, but the confirmation of these results is needed with genetic correlations over a larger dataset.

Spatio-temporal modeling of high-throughput multi-spectral aerial images improves agronomic trait genomic prediction in hybrid maize.

Design randomizations and spatial corrections have increased understanding of genotypic, spatial, and residual effects in field experiments, but precisely measuring spatial heterogeneity in the field remains a challenge. To this end, our study evaluated approaches to improve spatial modeling using high-throughput phenotypes (HTP) via unoccupied aerial vehicle (UAV) imagery. The normalized difference vegetation index (NDVI) was measured by a multi-spectral MicaSense camera and processed using ImageBreed. Contrasting to baseline agronomic trait spatial correction and a baseline multi-trait model, a two-stage approach was proposed. Using longitudinal NDVI data, plot level permanent environment (PE) effects estimated spatial patterns in the field throughout the growing season. NDVI PE were separated from additive genetic effects using two-dimensional spline (2DSpl), separable autoregressive (AR1) models, or random regression models (RR). The PE were leveraged within agronomic trait genomic best linear unbiased prediction (GBLUP) either modeling an empirical covariance for random effects, or by modeling fixed effects as an average of PE across time or split among three growth phases. Modeling approaches were tested using simulation data and Genomes-to-Fields (G2F) hybrid maize (Zea mays L.) field experiments in 2015, 2017, 2019, and 2020 for grain yield, grain moisture, and ear height. The two-stage approach improved heritability, model fit, and genotypic effect estimation compared to baseline models. Electrical conductance and elevation from a 2019 soil survey significantly improved model fit, while 2DSpl PE were most strongly correlated with the soil parameters. Simulation of field effects demonstrated improved specificity for RR models. In summary, the use of longitudinal NDVI measurements increased experimental accuracy and understanding of field spatio-temporal heterogeneity.

Comparison between multiple-trait and random regression models for genetic evaluation of weight traits in Australian meat sheep.

Random regression (RR) models are recommended as an alternative to multiple-trait (MT) models for better capturing the variance-covariance structure over a trajectory and hence more accurate genetic evaluation of traits that are repeatedly measured and genetically change gradually over time. However, a limited number of studies have been done to empirically compare RR over a MT model to determine how much extra benefit could be achieved from one method over another. We compared the prediction accuracy of RR and MT models for growth traits of Australian meat sheep measured from 60 to 525 d, using 102,579 weight records from 24,872 animals. Variance components and estimated breeding values (EBVs) estimated at specific ages were compared and validated with forward prediction. The accuracy of EBVs obtained from the MT model was 0.58, 0.51, 0.54, and 0.56 for weaning, postweaning, yearling, and hogget weight stages, respectively. RR model produced accuracy estimates of 0.56, 0.51, 0.54, and 0.54 for equivalent weight stages. Regression of adjusted phenotype on EBVs was very similar between the MT and the RR models (P > 0.05). Although the RR model did not significantly increase the accuracy of predicting future progeny performance, there are other benefits of the model such as no limit to the number of records per animal, estimation of EBVs for early and late growth, no need for age correction. Therefore, RR can be considered a more flexible method for the genetic evaluation of Australian sheep for early and late growth, and no need for age correction.

Currently, multiple-trait (MT) models are used in large-scale genetic evaluation of growth traits, where body weight traits are defined as separate traits at a finite number of fixed ages. Random regression (RR) models are expected to be superior since they can handle repeated measurements of weight and model these as a function of the actual age of measurement. These two models were compared in predicting breeding values for the body weight of Australian meat sheep. Phenotypic variation and estimated breeding values (EBVs) estimated at specific ages between 60 and 525 d with RR and MT models were compared and EBVs were validated in progeny data. The accuracy of EBVs in forecasting the performance of progeny was not statistically different between the two models. Other benefits of the RR model include the use of multiple records per animal, estimation of EBVs for early and late growth, with no need for age correction. Hence, RR models can be useful for the genetic evaluation of growth traits of sheep in Australia, but they do not necessarily predict breeding values at different ages more accurately than MT models.

B-spline polynomials models for analyzing growth patterns of Guzerat young bulls in field performance tests.

Objective: The aim of this study was to identify suitable polynomial regression for modeling the average growth trajectory and to estimate the relative development of the rib eye area, scrotal circumference, and morphometric measurements of Guzerat young bulls. Methods: A total of 45 recently weaned males, aged 325.8 ± 28.0 days and weighing 219.9 ± 38.05 kg, were evaluated. The animals were kept on Brachiaria brizantha pastures, received multiple supplementations, and were managed under uniform conditions for 294 days, with evaluations conducted every 56 days. The average growth trajectory was adjusted using ordinary polynomials, Legendre polynomials, and quadratic B-splines. The coefficient of determination, mean absolute deviation, mean square error, the value of the restricted likelihood function, Akaike information criteria, and consistent Akaike information criteria were applied to assess the quality of the fits. For the study of allometric growth, the power model was applied. Results: Ordinary polynomial and Legendre polynomial models of the fifth order provided the best fits. B-splines yielded the best fits in comparing models with the same number of parameters. Based on the -2RLL, AIC and CAIC criteria, the B-splines model with six intervals described the growth trajectory of evaluated animals more smoothly and consistently. In the study of allometric growth, the evaluated traits exhibited negative heterogeneity (b<1) relative to the animals' weight (p<0.01), indicating the precocity of Guzerat cattle for weight gain on pasture. Conclusion: Complementary studies of growth trajectory and allometry can help find interest when an animal's weight changes and thus assist in decision-making regarding management practices, nutritional requirements, and genetic selection strategies to optimize growth and animal performance.

Development of an index that decreases birth weight and promotes postnatal growth yet minimizes selection intensity in beef cattle.

Objective: The main purpose of our current study was to improve the growth curve of meat animals decreasing the birth weight yet achieves a finishing weight that is the same as that before selection but at younger age. Methods: Random regression model was developed to derive various selection indices to achieve desired gains in body weight at target time points throughout the fattening process. We considered absolute and proportional gains at specific ages (in weeks) and for various stages (i.e., early, middle, late) during the fattening process. Results: The point gain index was particularly easy to use because breeders can assign a specific age (in weeks) as a time point and model either the actual weight gain desired or a scaled percentage gain in body weight. Conclusion: The point gain index we developed can achieve the desired weight gain at any given postnatal week of the growing process and is an easy-to-use and practical option for improving the growth curve.

Modelling the lactation curve in Alpine × Beetal crossbred dairy goats using random regression models fitted with Legendre polynomial and B-spline functions.

The current study sought to genetically assess the lactation curve of Alpine × Beetal crossbred goats through the application of random regression models (RRM). The objective was to estimate genetic parameters of the first lactation test-day milk yield (TDMY) for devising a practical breeding strategy within the nucleus breeding programme. In order to model variations in lactation curves, 25,998 TDMY records were used in this study. For the purpose of estimating genetic parameters, orthogonal Legendre polynomials (LEG) and B-splines (BS) were examined in order to generate suitable and parsimonious models. A single-trait RRM technique was used for the analysis. The average first lactation TDMY was 1.22 ± 0.03 kg and peak yield (1.35 ± 0.02 kg) was achieved around the 7th test day (TD). The present investigation has demonstrated the superiority of the B-spline model for the genetic evaluation of Alpine × Beetal dairy goats. The optimal random regression model was identified as a quadratic B-spline function, characterized by six knots to represent the central trend. This model effectively captured the patterns of additive genetic influences, animal-specific permanent environmental effects (c2 ) and 22 distinct classes of (heterogeneous) residual variance. Additive variances and heritability (h2 ) estimates were lower in the early lactation, however, moderate across most parts of the lactation studied, ranging from 0.09 ± 0.04 to 0.33 ± 0.06. The moderate heritability estimates indicate the potential for selection using favourable combinations of test days throughout the lactation period. It was also observed that a high proportion of total variance was attributed to the animal's permanent environment. Positive genetic correlations were observed for adjacent TDMY values, while the correlations became less pronounced for more distant TDMY values. Considering better fitting of the lactation curve, the use of B-spline functions for genetic evaluation of Alpine × Beetal goats using RRM is recommended.

Longitudinal genetic dynamics of weaning index and implications for cow-calf production efficiency.

In beef cattle, the selection for higher weights at young ages has been questioned with the argument that this criterion may increase the adult weight of cows, resulting in higher costs. Therefore, selection criteria should be employed to increase weights at young ages with minimal impact on the adult weight of cows. Additionally, the relationship between measures of cow production efficiency and other well-established selection criteria in breeding programs remains poorly understood. The objective of this study was to longitudinally evaluate the relationship between the weaning index (WIndex) as a measure of efficiency and growth traits of the cows. Possible changes over time in WIndex due to selection applied for yearling weight (YW) were also investigated. The WIndex was proposed to maximize genetic response in the weaning weight of the calf while maintaining genetic gain in BW of the cow at zero. A random regression model was adopted to estimate correlations between WIndex, BW, hip height (HH), and body condition score (BCS) using records of Nelore cows from three lines. Genetic trends were calculated for the control line (NeC) and lines selected for greater YW (NeS and NeT). The age of 3 years was the most critical for the weaning efficiency of the cows. At this stage, young cows are still growing and wean lighter calves than their adult counterparts. The genetic correlation estimates between WIndex and BW (-0.58 to 0.04), HH (-0.05 to -0.34), and BCS (-0.51 to -0.17) were close to zero or negative. BW and HH were strongly correlated genetically across all ages (0.73-0.76). In general, HH exhibited a weak and negative genetic relationship with BCS. The genetic correlation between BW and BCS was stronger for advanced ages (0.45-0.68). In lines selected for YW, important increases in WIndex were observed. However, NeS has been selected since the 1980s until the present for YW, and thus, it showed a more pronounced trend of increasing BW and, consequently, a more modest trend of increasing WIndex compared to NeT. In contrast, WIndex exhibited a trend close to zero for NeC. In this context, monitoring HH and BCS can be useful to avoid losses in the weaning efficiency of cows. Furthermore, we suggest that one way to mitigate efficiency losses in calf production could involve stabilizing the BW of cows and increasing the weaning weight of calves using the WIndex.

Modeling daily yields of milk, fat, protein, and lactose of New Zealand dairy goats undergoing standard and extended lactations.

This study aimed to assess the milk production data for New Zealand dairy goats in either a standard lactation (SL; ≤305 d in milk [DIM]) or extended lactation (EL; >305 and ≤670 DIM) using a random regression (RR) with third- and fifth-order Legendre polynomials, respectively. Persistency of EL was defined as (B/A) × 100, where A was the accumulated yield from d 1 to 305, and B was the accumulated yield from d 366 to 670. On average, goats in SL produced 1,183 kg of milk, 37 kg of fat, 37 kg of protein, and 54 kg of lactose. The average production of milk, fat, protein, and lactose in EL were 2,473 kg, 78 kg, 79 kg, and 112 kg, respectively. The average persistences for milk, fat, protein, and lactose yields during EL were 98%, 98%, 102%, and 96%, respectively. The relative prediction errors were close to 10% and the concordance correlation coefficients >0.92, indicating that the RR model with Legendre polynomials is adequate for modeling lactation curves for both SL and EL. Total yields and persistency were analyzed with a mixed model that included the fixed effects (year, month of kidding, parity, and proportion of Saanen) as covariates and the random effects of animal and residual errors. Effects of year, month of kidding, and parity were significant on the total yields of milk, fat, protein, and lactose for both SL and EL. The total milk yield of first-parity goats with SL was 946 kg and the total milk yield of second-parity goats with SL was 1,284 kg, making a total of 2,230 kg over 2 years. The total milk yield of a first-parity goat with EL was 2,140 kg. Thus, on average, a goat with SL for the first and second parity produced 90 kg more milk than a first-parity goat subjected to EL. However, a second-parity goat subjected to EL produced 43 kg more milk (2,639 kg) than a goat with SL following the second and third parity (1,284 kg + 1,312 kg). These data, along with the various other benefits of EL (e.g., fewer offspring born and reduced risk of mastitis, lameness, and metabolic problems in early lactation), indicate that EL as a management strategy holds the potential to improve dairy goat longevity and lifetime efficiency without compromising milk production.

Candidate genes for longitudinal traits under sequential sampling in beef cattle.

Both the measurement age of a longitudinal trait and the common pre-sampling procedures used in beef cattle herds may affect the identification of a functional candidate gene (FCG) that is potentially associated with a trait. To identify the FCG that takes part in the genetic control of body weight at five different ages in a beef cattle population with and without sequential sampling, the animals were weighed at different measurement events, around 330, 385, 440, 495 and 550 days old. Genetic parameters were estimated for body weight at each age using a single trait (STM) and a random regression model (RRM). In addition, two different databases were used to estimate the genetic parameters: the first (DB100) was formed by all animals that were weighed in the five measurement events, and the second (DB70) has records of the same population, considering that 70% of the heaviest animals were selected after each measurement event. For DB100, genome-wide association studies (GWAS) were performed with 21,667 SNP markers to identify genomic windows that explained at least 1% of the genetic variance. Additionally, prioritization analyses were performed and FCGs were selected. We associated seven different FCGs with body weight at different ages. Among them, the gene DUSP10 was suggested as FCG in all five ages evaluated. Genetic parameters estimated for body weight using DB100 were similar when STM and RRM were applied. However, when DB70 was used as phenotypic data, there were differences between the two models. When the STM was applied, there were differences between the genetic parameters estimated for body weight when DB100 or DB70 were used as sources of phenotypes, but not for the estimates obtained with RRM. The importance of each gene for animal growth can change at different ages, and different genes may be more relevant to body weight at each different growth stage for beef cattle. Besides, sequential sampling can affect the GWAS results of a longitudinal trait. The age of the animal when a longitudinal trait is measured and pre-sampling can also contribute to inconsistencies in GWAS results for body weight in beef cattle, depending on the time when that data were collected, and consequently on the identification of FCG between studies, even when models that consider a covariance structure are used.

The marker effects of a single-step random regression model for 4 test-day traits in German Holsteins.

The single-step genomic model has become the golden standard for routine evaluation in livestock species, such as Holstein dairy cattle. The single-step genomic model with direct estimation of marker effects has been proven to be efficient in accurately accounting for millions of genotype records. For diverse applications including frequent genomic evaluation updates on a weekly basis, estimates of the marker effects from the single-step evaluations play a central role in genomic prediction. In this study we focused on exploring the marker effect estimates from the single-step evaluation. Phenotypic, genotypic, and pedigree data were taken from the official evaluation for German dairy breeds in April 2021. A multilactation random regression test-day model was applied to more than 242 million test-day records separately for 4 traits: milk, fat, and protein yields, and somatic cell scores (SCS). Approximately one million genotyped Holstein animals were considered in the single-step genomic evaluations including ∼21 million animals in pedigree. Deregressed multiple across-country breeding values of Holstein bulls having daughters outside Germany were integrated into the national test-day data to increase the reliability of genomic breeding values. To assess the stability and bias of the marker effects of the single-step model, test-day records of the last 4 yr were deleted, and the integrated bulls born in the last 4 yr were truncated from the complete phenotypic dataset. Estimates of the marker effects were shown to be highly correlated, with correlations ∼0.9, between the full and truncated evaluations. Regression slope values of the marker-effect estimates from the full on the truncated evaluations were all close to their expected value, being ∼1.03. Calculated using random regression coefficients of the marker effect estimates, drastically different shapes of the genetic lactation curve were seen for 2 markers on chromosome 14 for the 4 test-day traits. The contribution of individual chromosomes to the total additive genetic variances seemed to follow the polygenic inheritance mode for protein yield and SCS. However, chromosome 14 was found to make an exceptionally large contribution to the total additive genetic variance for milk and fat yields because of markers near the major gene DGAT1. For the first lactation test-day traits, we obtained ∼0 correlations of chromosomal direct genomic values between any pair of the chromosomes; no spurious correlations were found in our analysis, thanks to the large reference population. For trait milk yield, chromosomal direct genomic values appeared to have a large variation in the between-lactation correlations among the chromosomes, especially between first and second or third lactations. The optimal features of the random regression test-day model and the single-step marker model allowed us to track the differences in the shapes of genetic lactation curves down to the individual markers. Furthermore, the single-step random regression test-day model enabled us to better understand the inheritance mode of the yield traits and SCS (e.g., variable chromosomal contributions to the total additive genetic variance and to the genetic correlations between lactations).

Impact of heat stress on genetic evaluation of oocyte and embryo production in Gir dairy cattle.

Identifying and selecting genotypes tolerant to heat stress might improve reproductive traits in dairy cattle, including oocyte and embryo production. The temperature-humidity index (THI) was used, via random regression models, to investigate the impact of heat stress on genetic parameters and breeding values of oocyte and embryo production in Gir dairy cattle. We evaluated records of total oocytes (TO), viable oocytes (VO), cleaved embryos (CE), and viable embryos (VE) from dairy Gir donors. Twenty-four models were tested, considering age at ovum pick-up (AOPU) and THI means as a regressor in the genetic evaluation. We computed THI in eight periods, from 0 to 112 days before ovum pick-up, which were adjusted by different orders of Legendre polynomials (second, third, and fourth). The best-fit model according to Akaike's information criterion (AIC) and Model Posterior Probabilities (MPP) considered Legendre polynomials of third order and THI means of 112 days for TO, fourth order and 56 days for VO, second order and 28 days for CE, and second order and 42 days for VE, respectively. The heritability (h2) estimates across AOPU and THI scales ranged from 0.34 to 0.62 for TO, 0.31 to 0.58 for VO, 0.26 to 0.39 for CE, and 0.15 to 0.26 for VE, respectively. The fraction of the phenotypic variance explained by the permanent environment in different AOPU and THI scales ranged from 0.03 to 0.25 for TO, 0.05 to 0.26 for VO, 0.09 to 0.36 for CE, and 0.15 to 0.27 for VE, respectively. Spearman's rank correlation between the estimated breeding values in different AOPU and THI scale from the top 5% sires and females ranged from 0.18 to 0.90 for TO, 0.31 to 0.95 for VO, 0.14 to 0.85 for CE, and 0.47 to 0.94 for VE, respectively. The h2 estimates for all evaluated traits varied from moderate to high magnitude across AOPU and THI scales, indicating that genetic selection can result in rapid genetic progress for the evaluated traits. There was a reranking among the best animals in different AOPU and THI. It is possible to select dairy Gir cattle tolerant to heat stress to improve oocyte and embryo production.

Differential sensitivity of climate conditions on birth weight genetic values in mice divergently selected for birth weight residual variance.

After 32 generations of a divergent selection experiment for residual variance of birth weight in mice, two divergent lines were thus obtained: the heterogeneous line (H-line) and the homogeneous line (L-line). Throughout the generations, differences were observed between the two lines in traits such as litter size, survival at weaning, and birth weight variability caused by unidentified environmental conditions. The L-line exhibited advantages in terms of higher survival rates, larger litter sizes, and less sensitivity to changes in food intake. The study is an examination of the effects of climate as an environmental factor on the performance of these animals. Climate factors including maximum, minimum, and mean temperature (T), humidity (H), and TH index; at three stages (the fecundation, a week before the parturition and the parturition), were linked to a birth weight dataset consisting of 22,614 records distributed as follows: 8,853 corresponding to the H-line, 12,649 to the L-line, and 1,112 to the initial population. Out of the 27 analyzed climatic variables, the maximum temperature 1 wk before parturition (MXTW) was identified as the most influential when comparing heteroscedastic models with the deviance information criterion. The order of Legendre polynomial to apply in the following random regression model was tested by a cross-validation using homoscedastic models. Finally, MXTW was compared on how it affected the two divergent lines by analyzing predicted breeding values (PBV) obtained from a random regression heteroscedastic model. The mean PBV of the H-line in the first generation showed a range of 0.070 g with a negative slope, which was 35 times higher than the range obtained for the L-line, which varied within 0.002 g. In the last generation of selection, the H-line exhibited greater instability of PBV across temperatures, with a difference of 0.101 g between the maximum and minimum mean PBV, compared to 0.017 g for the L-line. The standard deviations of the slopes in the H-line were more dispersed than in the L-line. Unlike the H-line, the L-line had slopes that were not significantly different from 0 throughout the generations of selection, indicating greater stability in response to MXTW variations. The H-line exhibited a higher sensitivity to changes in MXTW, particularly in birth weight, with the L-line being more stable. The selection for uniformity of birth weight could lead to less sensitive animals under environmental changes.

Two mice lines obtained by divergent selection for birth weight residual variance were used to determine whether environmental factors could differently affect the homogeneous and heterogeneous lines. The maximum temperature 1 wk before parturition (MXTW) had the higher impact on the birth weight of the animals. A random regression model showed the individual trajectory of birth weight throughout the changes in MXTW. It was evident that the homogeneous line is less susceptible to changes in climate. This result, therefore, supports the hypothesis that the selection for homogeneity in production animals is more advantageous. More robust animals are obtained that can better cope with changes in climate without compromising their productive traits.

Novel genetic parameters for genetic residual feed intake in dairy cattle using time series data from multiple parities and countries in North America and Europe.

Residual feed intake is viewed as an important trait in breeding programs that could be used to enhance genetic progress in feed efficiency. In particular, improving feed efficiency could improve both economic and environmental sustainability in the dairy cattle industry. However, data remain sparse, limiting the development of reliable genomic evaluations across lactation and parity for residual feed intake. Here, we estimated novel genetic parameters for genetic residual feed intake (gRFI) across the first, second, and third parity, using a random regression model. Research data on the measured feed intake, milk production, and body weight of 7,379 cows (271,080 records) from 6 countries in 2 continents were shared through the Horizon 2020 project Genomic Management Tools to Optimise Resilience and Efficiency, and the Resilient Dairy Genome Project. The countries included Canada (1,053 cows with 47,130 weekly records), Denmark (1,045 cows with 72,760 weekly records), France (329 cows with 16,888 weekly records), Germany (938 cows with 32,614 weekly records), the Netherlands (2,051 cows with 57,830 weekly records), and United States (1,963 cows with 43,858 weekly records). Each trait had variance components estimated from first to third parity, using a random regression model across countries. Genetic residual feed intake was found to be heritable in all 3 parities, with first parity being predominant (range: 22-34%). Genetic residual feed intake was highly correlated across parities for mid- to late lactation; however, genetic correlation across parities was lower during early lactation, especially when comparing first and third parity. We estimated a genetic correlation of 0.77 ± 0.37 between North America and Europe for dry matter intake at first parity. Published literature on genetic correlations between high input countries/continents for dry matter intake support a high genetic correlation for dry matter intake. In conclusion, our results demonstrate the feasibility of estimating variance components for gRFI across parities, and the value of sharing data on scarce phenotypes across countries. These results can potentially be implemented in genetic evaluations for gRFI in dairy cattle.

Analysis of Early Growth of Piglets from Hyperprolific Sows Using Random Regression Coefficient.

Management of hyperprolific sows is challenging when it comes to controlling birth weight variability and piglet survival in large litters. The growth of low birth weight piglets can be compromised and have a negative impact on production efficiency. The objective of the study was to apply a random regression coefficient model to estimate the main effects of the growth of piglets of highly prolific sows. The dataset contained growth data for 360 piglets from 25 Pen Ar Lan Naima sows. In addition to routine procedures after farrowing, piglets were weighed five times: on day 1 after farrowing, on day 14 of life, at weaning on day 28, on day 30 of nursery period, and at the end of the nursery period when piglets were 83 days old. Data were treated as longitudinal, with body weight as the dependent variable. Fitting age as a quadratic regression within piglets in the random part of the model helped to determine the significant effect of birth weight, litter size, and parity on the growth of the piglets. Since the piglets from large litters often have non-uniform birth weights and this can affect further growth, the use of a random regression coefficient model is practical for analysing the growth of such piglets due to the ability to describe the individual growth pattern of every individual.

Estimating breeding values for feed efficiency in dairy cattle by regression on expected feed intake.

The efficiency with which a dairy cow utilises feed for the various physiological and metabolic processes can be evaluated by metrics that contrast realised feed intake with expected feed intake. In this study, we presented a new metric - regression on expected feed intake (ReFI). This metric is based on the idea of regressing DM intake (DMI) on expected DMI using a random regression model, where energy requirement formulations are applied for the calculation of expected DMI covariables. We compared this new metric with the metrics residual feed intake (RFI) and genetic residual feed intake (gRFI), by applying them on 18 581 feed efficiency records from 654 primiparous Nordic Red dairy cows. We estimated variance components for the three metrics and their respective genetic correlations with intake and production traits. In addition, we examined the phenotypes of superior cows. With ReFI, we estimated for feed efficiency a higher genetic variation (4.7%) and heritability (0.23) compared to applying RFI or gRFI. The ReFI metric was genetically uncorrelated with DMI and negatively correlated within energy-corrected milk (ECM), whereas the RFI metric was genetically positively correlated with DMI and metabolic BW. The gRFI metric was genetically positively correlated with DMI and uncorrelated with energy sink traits. Overall, the estimated SE were large. The ReFI metric resulted in a different ranking of cows compared to those based on RFI or gRFI and was superior in selecting the most efficient animals. When the selection was based on ReFI breeding values, then the 10% most efficient cows produced 12.3% more ECM per unit metabolisable energy intake, whereas the corresponding values were only 4.3 or 5.9% when using RFI or gRFI breeding values, respectively. Based on ReFI, superior cows had also higher milk production, whereas based on RFI or gRFI milk production either decreased or was unaffected, respectively. The superiority of the ReFI metric in selecting efficient cows was due to a better modelling of the expected feed intake. The ReFI metric simplified modelling of feed utilisation efficiency in dairy cattle and resulted in breeding values that are equal to percentages of feed saved.

Random regression for modeling soybean plant response to irrigation changes using time-series multispectral data.

Plant response to drought is an important yield-related trait under abiotic stress, but the method for measuring and modeling plant responses in a time series has not been fully established. The objective of this study was to develop a method to measure and model plant response to irrigation changes using time-series multispectral (MS) data. We evaluated 178 soybean (Glycine max (L.) Merr.) accessions under three irrigation treatments at the Arid Land Research Center, Tottori University, Japan in 2019, 2020 and 2021. The irrigation treatments included W5: watering for 5 d followed by no watering 5 d, W10: watering for 10 d followed by no watering 10 d, D10: no watering for 10 d followed by watering 10 d, and D: no watering. To capture the plant responses to irrigation changes, time-series MS data were collected by unmanned aerial vehicle during the irrigation/non-irrigation switch of each irrigation treatment. We built a random regression model (RRM) for each of combination of treatment by year using the time-series MS data. To test the accuracy of the information captured by RRM, we evaluated the coefficient of variation (CV) of fresh shoot weight of all accessions under a total of nine different drought conditions as an indicator of plant's stability under drought stresses. We built a genomic prediction model (MTRRM model) using the genetic random regression coefficients of RRM as secondary traits and evaluated the accuracy of each model for predicting CV. In 2020 and 2021,the mean prediction accuracies of MTRRM models built in the changing irrigation treatments (r = 0.44 and 0.49, respectively) were higher than that in the continuous drought treatment (r = 0.34 and 0.44, respectively) in the same year. When the CV was predicted using the MTRRM model across 2020 and 2021 in the changing irrigation treatment, the mean prediction accuracy (r = 0.46) was 42% higher than that of the simple genomic prediction model (r =0.32). The results suggest that this RRM method using the time-series MS data can effectively capture the genetic variation of plant response to drought.

Antagonistic Effects of Assortative Mating on the Evolution of Phenotypic Plasticity along Environmental Gradients.

AbstractPrevious theory has shown that assortative mating for plastic traits can maintain genetic divergence across environmental gradients despite high gene flow. Yet these models did not examine how assortative mating affects the evolution of plasticity. We here describe patterns of genetic variation across elevation for plasticity in a trait under assortative mating, using multiple-year observations of budburst date in a common garden of sessile oaks. Despite high gene flow, we found significant spatial genetic divergence for the intercept, but not for the slope, of reaction norms to temperature. We then used individual-based simulations, where both the slope and the intercept of the reaction norm evolve, to examine how assortative mating affects the evolution of plasticity, varying the intensity and distance of gene flow. Our model predicts the evolution of either suboptimal plasticity (reaction norms with a slope shallower than optimal) or hyperplasticity (slopes steeper than optimal) in the presence of assortative mating when optimal plasticity would evolve under random mating. Furthermore, a cogradient pattern of genetic divergence for the intercept of the reaction norm (where plastic and genetic effects are in the same direction) always evolves in simulations with assortative mating, consistent with our observations in the studied oak populations.

Relationship between milk urea content and important milk traits in Holstein cattle.

Breeding cattle with low nitrogen emissions has been proposed as a countermeasure against eutrophication due to dairy production. Milk urea content (MU) could potentially serve as a new readily measured indicator trait for nitrogen emissions by cows. Therefore, we estimated genetic parameters related to MU and its relationship with other milk traits. We analysed 4 178 735 milk samples collected between January 2008 and June 2019 from 261 866 German Holstein dairy cows during their first, second, and third lactations. Restricted maximum likelihood estimation was conducted using univariate and bivariate random regression sire models in WOMBAT. We obtained moderate average daily heritability estimates for the daily MU of 0.24 in first lactation cows, 0.23 in second lactation cows, and 0.21 in third lactation cows with average daily genetic SDs of 25.16 mg/kg, 24.93 mg/kg, and 23.75 mg/kg, respectively. Averaged over days in milk, the repeatability estimates were low at 0.41 in first, second, and third lactation cows. A strong positive genetic correlation was found between MU and milk urea yield (MUY; 0.72 on average). In addition, 305-day heritabilities were estimated as 0.50, 0.52, and 0.50 in first, second, and third lactation cows, respectively, with genetic correlations of 0.94 or higher for MU in different lactations. By contrast, the averaged estimates of the genetic correlations between MU and other milk traits were low (-0.07 to 0.15). Moderate heritability estimates clearly allow the possible selection for MU, and the near-zero estimates of genetic correlations indicate no risk of undesired correlated selection responses in other milk traits. However, a relationship still needs to be established between MU as an indicator trait and the target trait, defined as total individual nitrogen emissions.